Extracting the Q electric quadrupole moment from lattice QCD data 
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The Q,~ has an extremely long lifetime, and is the most stable of the baryons with spin 3/2. 
Therefore the fl~ magnetic moment is very accurately known. Nevertheless, its electric quadrupole 
moment was never measured, although estimates exist in different formalisms. In principle, lattice 
QCD simulations provide at present the most appropriate way to estimate the Q~ form factors, as 
function of the square of the transferred four- momentum, , since it describes baryon systems at 
the physical mass for the strange quark. However, lattice QCD form factors, and in particular Ge2, 
are determined at finite only, and the extraction of the electric quadrupole moment, Q^- = 
Ge2{0) 2Mii ' involves an extrapolation of the numerical lattice results. In this work we reproduce 
the lattice QCD data with a covariant spectator quark model for Q,^ which includes a mixture of S 
and two D states for the relative quark-diquark motion. Once the model is calibrated, it is used to 
determine Qn-- Our prediction is Qn- = (0.96 ± 0.02) x 10"^ efm^ [Ge2(0) = 0.680 ± 0.012]. 
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I. INTRODUCTION 

The prediction of the electromagnetic structure of 
baryons and mesons is an important challenge for quark 
models and, when compared to the available experimen- 
tal results, it provides a test on the relevant hadronic 
degrees of freedom. For baryons with spin 1/2 (as in 
the baryon octet) the charge and the magnetic moment 
are the only multipole moments to be defined, while 
for baryons with spin 3/2 (as in the baryon decuplet) 
also the electric quadrupole moment exists. However, 
at present, there is no experimental measurement of the 
electric quadrupole moment for any of the baryons, al- 
though there are several model predictions for the A and 
the ri~, and other decuplet particles. For a summary of 
the A results see Refs. [E 111- As for the ri~ there are 
predictions based in quark models [3i-t9|L chiral pertur- 
bation theo ry j iol. large- A^c limit |13 - [T^ and other 
formalisms |15l4l8l |. 

Within all members of the baryon decuplet, the il^ is 
specially interesting. As it is composed solely by strange 
quarks in the valence sector, it can decay only by weak 
interaction and therefore its lifetime is extremely longer 
than the one of the other baryons. For this reason, exper- 
imentally, the r2~ properties are easier to be determined 
than the ones of any other member of the decuplet. A 
good illustration of this is the accuracy of the Vt~ mag- 
netic moment — —(2.019 ± 0.053)^Ar p^l - [2^ . where 
fiN the nuclear magneton. It is also expected that 
the Qo - q uadrupole moment will be measured in a near 
future dJii]. 

Additionally, the valence quark content of the is 
restricted to strange quarks, with a mass considerably 
larger than the light u and d quark masses. As now it is 
already possible to perform lattice QCD at the physical 
strange quark mass, the f2~ magnetic moment [27H29|, 



and more recently also its electric charge Geo, magnetic 
dipole Gmii electric quadrupole Ge2 and magnetic oc- 
tupole Gm3 [M, HH were calculated within lattice QCD. 
Another important issue is that in sea quark effects for 
the only at most one single light quark participates, 
and therefore the pion has no role in this case. As in 
chiral perturbation theory loops involving mesons heav- 
ier than the pion are suppressed, the 57" becomes then a 
special case where meson cloud corrections to the valence 
quark core are expected to be small. A consequence of 
the smallness of the meson cloud effects is that lattice 
QCD simulations, quenched or unquenched, should be a 
good approximation to form factors at the physical 
point. Therefore in this work we take the lattice QCD 
simulations as good representations of the physical re- 
sults, without any extrapolation of the lattice data to 
the physical pion mass. 

The main limitation in obtaining the f2~ electromag- 
netic form factors in lattice QCD simulations comes from 
these ones being restricted, for practical reasons, to fi- 
nite non-zero values of Q^, while the determination of 
the quadrupole moment, for instance, is proportional to 
Ge2{Q)- An extrapolation in the momentum transfer 
squared Q^, down to = is then required, and one 
has inevitably to resort to an analytical form to do it. 
It is at this point that it is reasonable to expect that a 
quark model is useful, in particular for systems as the 
f2~, without light valence quarks and where meson loop 
corrections are expected to be small. Since the covariant 
spectator quark model was tested already for spin 1/2 
baryons [Sa - fsst . spin 3/2 baryons [U, 0, 0, [H, 1131, in- 
cluding strange quarks, and electro mag netic transitions 
between different baryon states [H, |37H4^ , it is a good 
candidate to be used not only to interpolate between lat- 
tice QCD data in Q^, but also to extrapolate the form 
factor data to = 0. 



2 



The covariant spectator quark model was applied in 
the past to estimate successfully the leading order form 
factors of the fl^ {Geo and Gmi) neglecting D-state ad- 
mixtures [131 • Here, we extend the formalism to the 
case where the D-state admixture coefficients are non- 
zero. Because these states are included, we have now 
contributions for the Ge2 and Gms form factors, and 
therefore we can use the model to extract the n,~ electric 
quadrupole moment from lattice QCD data. By adjust- 
ing some parameters associated with the fl^ wavefunc- 
tion to the lattice data, our procedure has the advan- 
tage of incorporating into the phenomenological model 
the fundamental theory of the strong interaction, in its 
discrete version (lattice QCD). The information on the 
wavefunction parameters allows us then to calculate all 
the electromagnetic form factors Geo, Gmi, Ge2 and 
Gm3 as functions of . In particular the model can be 
used to determine the electric quadrupole form factor at 
the = point. To constrain the model we use the 
unquenched lattice QCD data from Ref. [s^ for at 
the physical $7" mass. Although the existing lattice data 
is unquenched, the valence quarks are expected to play 
the main role and meson dressing to be small, as men- 
tioned before. This is why the adjustment of the quark 
model to the lattice QCD data is meaningful. To better 
constrain the parameters of the quark model, we use also 
the single datapoint for Gm3 from Ref. [3l|, in addition 
to the lattice data from Ref. [30| . 

We start by using the spectator formalism to repre- 
sent the J7~ wavefunction, in a similar way to the one 
used before for the A Although both the A and 

the VL~ have the same spin structure, they differ in their 
flavor content. Therefore we begin with the SU(3) gen- 
eralization of the spectator quark model for the overall 
study of the baryon decuplet [13]. Also, the quark mo- 
mentum distributions for the two baryons are different, 
and our calculation at the end quantifies this difference. 
As for the electromagnetic current associated with the 
interaction of the photon with the strange quark, which 
is another aspect where the calculation differs from the 
calculation for the A baryon, we use the current based 
on vector meson dominance from Ref. [20j . 

Under the assumption that the D-state components 
are small, then we take only the electromagnetic current 
matrix elements which are in first order in the admixture 
coefficients as done in Ref. [2] for the A. Finally, in the 
process of adjusting the electromagnetic form factor 
results to the lattice data, we determine the percentage 
of each D state present in the orbital quark-diquark part 
of the wavefunction of the baryon. Our calculation 
enables us at the end to narrow the uncertainty in the ex- 
traction of the value of the quadrupole magnetic moment 
of the 51^ from the lattice QCD data. 

This work is organized as follows: In Sect. |n]we give 
the formulas for the il^ electromagnetic form factors in 
first order of the admixture coefficients; In Sect. IIIII we 
parametrize the wavefunction and its momentum de- 
pendence. The results are presented in Sect. II VI and the 



conclusions and final remarks in Sect. [V] 

II. Q- FORM FACTORS 

We use here the covariant spectator quark model, 
where relativity is implemented consistently. Within this 
framework a baryon is described described as a off-shell 
quark and two noninteracting on-shell spectator quarks. 
Integrating over the on-mass-shell quarks degrees of free- 
dom one represents those quarks states as a single on- 
shell particle (or diquark) with an average mass m^) 
[20I |32| . With this reduction, the wavefunction associ- 
ated with the baryon states including the spin, angular 
momentum, coordinate space and flavor structure, can 
then be written as the direct product of the diquark and 
quark states properly symmetrized. 

The electromagnetic interaction with the n~ is, in 
relativistic impulse approximation, written as the sum 
over the terms in which the photon couples to each (off- 
shell) quark in turn with the other two (on-shell) quarks 
that compose the diquark. The electromagnetic struc- 
ture of the quarks is parametrized in terms of form fac- 
tors. One can re-arrange the contributions to the elec- 
tromagnetic current from the conveniently symmetrized 
wavefunctions in terms of the on-shell diquark states - 
quark pair (12), and the off-shell quark - quark 3 states. 
The final result for the current becomes then three times 
the current associated to the interaction with quark 3. 
See Ref. [13] for a detailed discussion. We write then 
the fl~ wavefunction as \E'a(P, fc) for total momentum 
P, as the combination of the diquark (on-shell) states, 
with momentum fc, and the quark 3 (off-shell) states. 
In our notation the indices for the diquark polarization 
e = 0, ± and the spin projection are omitted for sim- 
plicity, and the matrix element of the electromagnetic 
current between the initial and final states of momentum 
P_i_ and P_ , respectively, is written as 

J'' = / ^niP+,k)j!;^n{P-,k), (1) 

g J k 

where jj^ is the current operator for quark 3 and J^, is the 
covariant integral is defined as = / (27r)^2gc ' where 
Ed is the diquark on-shell energy. In Eq. ([T|) the interac- 
tions with all quarks are counted, without including the 
coupling with the diquark [20l . |32| . 

In this work the 51^ wavefunction is represented as a 
combination of an S state and two D states for the quark- 
diquark relative motion (39j : 

^niP, k)^N [*s(P k) + a^^DsiP, k) + 6*di(P, k)] . 

(2) 

In the previous equation a and h are the mixture coeffi- 
cients of the states: D3 (core spin 3/2) and Dl (core spin 
1/2) respectively, and N = l/Vl -f -|- 6^ ^ normal- 
ization constant. Each of the three wavefunction terms 
includes a scalar wavefunction, respectively ips, ipDs and 
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ipDi, which is a function of the baryon momentuni P 
and diquark momentum k, with a form and normahza- 
tion given in the next section. 

The quark current jj^ can be in general decomposed as 

m')-.niQ'h'+MQ')'-^, (3) 



order in the admixture coefficients a and b: 



2Mn ' 



where Mm is the nucleon mass, and 
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(4) 



with Ao = diag(f,l,0): A3 = diag(l,— 1,0) and A., = 
diag(0, 0, —2) are SU(3) flavor operators acting in the 
particular quark flavor state q = {uds)^ [20I, [3^. 

The quark form factors are normalized as /i±(0) = 1, 
/io(0) = 1, /2±(0) = K± and /2o(0) = k,. We model 
the electromagnetic structure of the quarks by means of 
a parametrization that is based on vector meson domi- 
nance. In particular the strange quark form factors are 
represented [l^l by 



/io = A+(l-A)- 



m5 



- Co 



/20 = { do-zr^ + (1 - do)- 



(5) 



(6) 



where is a mass of the (j) meson (system ss), Mh — 
2Mjv is an effective vector meson that simulates the short 
range structure and A is fixed by deep inelastic scattering, 
and corresponds to the quark number density (32j. The 
coefficients defining the current were fixed as cq = 4.427 
and do = —1-860 in the study of the dominant form fac- 
tors of the baryon decuplet ^] . As for the strange quark 
anomalous magnetic moment one uses Kg = 1.462 to re- 
produce the experimental value of /i^- [13] ■ The explicit 
expression for the remaining quark form factors are pre- 
sented in Refs. (Hil El- 

Following Refs. [1^ H^] the current can be written in 
terms of charge en and anomalous magnetic moment ko 
functions 





= N^gnXs 




(11) 


Gmi{Q^) 
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(12) 


Gb2(Q') 


= N^-gni-ia)^ 

T 




(13) 


Gm3,{Q^) 




T T 
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(14) 



where the overlap between the S-states is 
Is^ J^4's{P+,k)i;siP-,k), 



(15) 



and the overlap between the S and each one of the D- 
states is 

Idi = J^bCk+,q+)^Di{P+,k)^Js{P-,k) (16) 
Id3^ f bik+,q+)i;D3{P+,k)ijsiP-,k). (17) 

Jk 

The function 6(fc+, 9+) is defined in Ref. [s^ and includes 
the specific angular dependence of a D-state, with 
defined as the three-momentum in the frame where P+ — 
(A/0,0, 0,0): k+^k~^P+. 

From Eq. pT|) one obtains that the result for the charge 
form factor at = is G'_bo(0) = — 7V^, [note that 
N"^ = l/(l + a^ + 6^)], which differs from the exact result 
(—1), if a, 6 7^ 0. This deviation, although small if a and b 
are small, is a consequence of taking in the calculation of 
the current matrix elements only the terms in first order 
in these D state admixture coefficients a and b. Also 
the magnetic dipole form factor Gmi, from Eq. (|12p . at 
— deviates slightly from the experimental value 
(which gives the magnetic moment). Once the terms 
for the current matrix elements for the D- to D-state 
transitions are included, the exact results of both the 
charge and magnetic moment are recovered exactly. We 
will use this fact to estimate the error of our model, as 
explained in Sect. IIVI 



en 



^/lo(Q') 



(7) 

(8) 



In the previous equation, for convenience we use a nor- 
malization that differs from the one presented in (20| . 
by a factor jj^. This redefinition does not change the 
results. We define also 

go = eo - TKQ (9) 
fn = en + kn, (10) 

where T = QV(4Mj2^). 

Working the algebra for the current as in Refs. [H, 0! 
replacing Sa and /ca by and fco, one obtains in first 



III. MODEL FOR THE SCALAR 
WAVEFUNCTIONS 

To describe the momentum dependence of the scalar 
wavefunctions one assumes a certain form (as done al- 
ready for the A in Refs. [li,!!^]) 



^s{P,k)^ 
^Dz{P.k) -- 
^Di{P,k) = 



Ns 



moias + Xn)^ 



Ndi 



(18) 
(19) 
(20) 
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where 



Xr, 



Mnmo 



(21) 



In this way we introduce momentum scale parameters 
(as,aD3,C(Di) for each angular momentum or orbital 
state. In contrast to the A case one does not have to 
impose the orthogonality of the spin core S* = 1/2 D- 
state with an S-state with the same spin, as a spin 1/2 
S-state with three strange quarks is forbidden. There- 
fore, the expression for the state Dl differs from the one 
used for the A baryon [H, , and one may consider a 
simpler parameterization for the fl~ Dl state. 
The normalization conditions are given by 



|P^m(-P,fc)P = l, 



(22) 
(23) 
(24) 



where P = (Mq, 0,0,0) is the momentum of the ft^ in 
the rest frame. The parameters as,aD3,ctDi and the 
admixture coefficients a, b are the five adjustable param- 
eters of the model. 



IV. RESULTS 

In this section we present the results for the ft~ electro- 
magnetic form factor data obtained from the quark model 
described in the previous sections, calibrated by a fit to 
the available lattice data. The available f2~ form factor 
data as function of Q^, is restricted to the unquenched 
lattice simulations of the Geo, Gaii and Ge2 form fac- 
tors [30i]. Under these conditions, and given that the 
term from the S-state dominates in Geo and Gmi, the 
most important constraint for the D-states comes from 
Ge2 (to which only the interference term between S- and 
D3-state contributes), and a decisive constraint on the 
Dl state (through Gms) is not available yet, except for 
the quenched calculation of Ref. [U for = o.23 GeV^. 
That point is also considered in our fit. 

In principle one could assume that the Dl state would 
not be important for the structure — we evoke that 
it was checked earlier that the lattice data for the A 
baryon can be described with very small (« 1%) D state 
admixtures Still, this needs to be investigated, and 
therefore we check here whether the Dl state can be im- 
portant to Geo and Gmi data and helps to improve the 
overall description of the form factor data. This is why in 
this work we leave the Dl-state mixture free and use the 
lattice data (mainly Gmi) to fix that contribution. Fu- 
ture lattice QCD simulations for Gm^ may then confirm 
or contradict our model. 

The only limitation of our calculation is that the 
amount of D-states admixture is assumed to be small. 



since the calculation of the form factor proceeds by tak- 
ing only into account the first order terms in a and 6 
(neglecting transitions between D-states). That limita- 
tion will be quantified at the end, comparing the result 
of G£;o(0) with the fi" charge (—1), obtained when all 
states are include in the electromagnetic transition cur- 
rent. 



A. Lattice data 

We use the lattice QCD data from Alexandrou et 
al. and the Gmz datapoint from Ref. [3l|. All the 
lattice QCD simulations from Ref. (30i] are unquenched 
but two different methods are used: hybrid action and 
domain wall fermions (DWF). The single datapoint from 
Ref. [3l| is quenched. The data from Ref. Q goes al- 
most to 4 GeV2 dl, but only the data for < 1.6 GeV^ 
is relatively precise. The simulations were performed for 
one pion mass using the hybrid method (m^r — 0.353 
GeV) and three pion masses with the DWF method 
(m^ = 0.297, 0.330 and 0.355 GeV). 

The results of Geo and Gmi in the DWF case show a 
weak dependence in the pion mass, suggesting that me- 
son cloud effects {K and ??) may be negligible at least for 
those form factors as expected from a three strange va- 
lence quark system. The hybrid simulation for Geo and 
Gmi shows a systematic deviation from the results from 
DWF (slower falloff) for large with similar pion mass 
(to^ = 0.355 GeV). 

For the hybrid simulation one has results for Geo, Gmi 
and Ge2 for pion mass value of — 0.353 GeV. As 
for the DWF simulations, one has data for Gmi, with 
TOtt — 0.355 GeV; for Geo and Gmi with m^r = 0.330 
GeV; and finally for data for Geoi Gmi and also for Ge2 
with = 0.297 GeV. The difference between the two 
methods, particularly to high Q^, can be a lattice artifact 
(cutoff effect) [33|. More lattice QCD simulations with 
smaller lattice spacing and larger volumes are necessary 
to clarify the differences between the two methods ;iO . 

The form factors can also be obtained from a sim- 
ulation of the A^ in the SU(3) symmetry limit. This was 
done by Boinepalli et al. [31| for pion mass — 0.697 
GeV, leading to a slightly heavier mass for the Vl~ (1.732 
GeV to be compared with 1.672 GeV of the physical case 
p^). This simulation is quenched and the results are 
given only for one value of [Q^ = 0.23 GeV^]. Impor- 
tantly, that work provides the only existing clue for the 
behavior of the octupole magnetic moment form factor: 



Gm3(0.23 GeV^) = 1.25 ± 7.50. 



B. Calibrating the model (fit to lattice data) 

Ideally, to extrapolate the form factors from lattice 
QCD one should take the set of lattice QCD simulations 
performed as close as possible to the physical limit - the 
lowest pion mass (to^ = 0.297 GeV) considered in the 
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x'(GAfi) 


x'(Ge2) 


x'(Gm3) 


x'M) Gb2(0) 


5.28 


2.66 


0.203 


2.49 


2.58 0.674 



TABLE I: Results of the in fit to the lattice data. 
The results correspond to the values: a =0.034f , b = 0.2666, 
as = 0.f793, am = 0.5394 and am = 0.4674. 



DWF case. Unfortunately, the available DWF data for 
Ge2 is restricted to 7 data points, for that mass, below 2 
GeV^, and there is no data for < 0.4 GeV^, which is a 
severe limitation to extrapolate the results down to = 
0. To obtain a more accurate extrapolation, we took also 
the hybrid data (with 13 datapoints for Ge2 below 2 
GeV^). To consider simultaneously the hybrid and DWF 
data increases the statistics, which in principle leads to 
more accurate constraint of the model. However, the 
differences between the two methods, in particular the 
high accuracy of the Geo data, make a good fit difficult, 
particularly for the Geo and Gmi form factors. Since 
the difference between the two methods is amplified as 

increases (particularly for > I GeV^), to improve 
the quality of our fit we took only the Geo and Gmi 
data below < 1 GeV^. This is justified since the 
lattice simulations have larger numerical error for higher 

and also because we are focused in the extrapolation 
for the — point, which leads us to reduce the weight 
of high data. This procedure is frequently used in 
lattice QCD studies at low The inclusion of the 

region > 1 GeV^ in the fit will be possible once a 
more homogeneous set of data with sufficient statistics is 
provided. As for Ge2 we took the data for < 2 GeV^, 
since the errorbars are more significant and we want to 
keep the statistics as large as possible. 

The results for the of tti^ fit are presented in the 
Table m The large for the Geo data is a consequence 
of the great precision of the data and also of the different 
behavior of the hybrid and DWF data. The parameters 
of the fit are indicated in the caption of Table HI The 
mixture of D-states that we obtain is 0.11% (for D3) and 
6.6% (for Dl). There is therefore a significant mixture of 
the Dl state. 

In detail, the value for the momentum-range parameter 
as (0.1793) obtained in this fit (with S and D states) 
suggests that D states improve in fact the description of 
Geo and Gmi, when compared with the predictions from 
Ref. [20] , prior to the simulations of Ref. [30] and where 
as was smaller (0.1630). 

The results of our fit to the lattice form factors are 
presented in Fig. [TJ The upper (lower) results for Geo 
and Gmi {Ge2 and Gm^) are given by Eqs. (ITT|) -(ITil). 
with the factor iV^ = + + b'^) defined by the 
normalization of the 17^ wavefunction^, given by Eq. 



^ Note that in the present work, by taking A^^ = 1/(1 + a? + b^), 



As mentioned in Sect. |TI1 those equations are derived in 
first order in the D-state admixture parameters a and b, 
and differ from the final result because the transitions 
between D-states are neglected. We can calculate the 
effect of this approximation by evaluating the correction 
to Geo{Q^) needed to reproduce exactly the ft~ charge 
at (5^ =0. From Eq. ([fl]) the exact result for G_bo(0) 
is obtained by setting N'^ = 1, which then defines our 
lowest estimation of Geo- The band between the upper 
and the lower result measures the effect in the final result 
of the neglected D to D-state transitions. The width of 
the band is small because a and b turn out to be small 
from the fit to the lattice data. As for the uncertainty in 
the calculation of the form factors, Gmi, Ge2 and Gms 
we use the same procedure, for consistency. 

Our results presented in Fig. [TJ for the four form fac- 
tors, are consistent with the overall lattice data within 
one standard deviation, to the exception of the form fac- 
tor Gm3- In that case, although we overestimate the lat- 
tice datapoint [Slj , our result is still inside an uncertainty 
interval of 1.5 standard deviations. Note that in our for- 
malism the magnitude of Gais is a consequence of the 
relatively large Dl state admixture (6.6%), a prediction 
to be tested by future and more precise lattice QCD sim- 
ulations. The extracted value of the electric quadrupole 
form factor at = q is Ge2{0) = 0.680 ± 0.012. 

V. CONCLUSIONS 

The 51^ is the most stable baryon with spin 3/2. Yet, 
only some of its properties are known. As a spin 3/2 par- 
ticle with a long lifetime, it is the first candidate for the 
experimental determination of the electric quadrupole 
moment, since the nucleon (spin 1/2) has no quadrupole 
moment, and the lifetime of the A is much shorter. Sev- 
eral experimental methods have been proposed to mea- 
sure the still unknown quadrupole moment Qn-, which 
is a signature of distortion and is likely to be measured 
in a very near future. This makes its a-priori prediction 
so challenging. 

As the ri~ is essentially a three strange valence quark 
system it is possible nowadays to simulate the elec- 
tromagnetic coupling with the $7" in a discrete lattice 
for physical strange quark masses with light sea quarks 
(tott — 300 MeV) and to determine the n~ electromag- 
netic form factors. In the absence of experimental in- 
formation, lattice QCD provides then the more reliable 
method to unveil the n~ electromagnetic structure. 

However, to determine the electric quadrupole moment 
an extrapolation of the G e2 form factor down to = 



we are not modifying the normalization of the H ^ wavefunction 
relatively to Ref. [M . In previous calculations we had A'^'^ = 1 
because only the S state (normalized to 1) was considered. Once 
the D- states are added, A^^ is redefined accordingly to the D- 
state admixture coefficients to N'^ = 1/(1 + -1-6^). 
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FIG. 1: Best fit with a Dl and D3 mixture. Lattice QCD data from Alexandrou et al. [30ll . For Gmi we include also the experimental 
result Gmi(0) = -3.604 ± 0.096 from PDG fill (*). The open circles represent the result for = o.23 GeV^ from Boinepalli et al. [sj]. 



must be done, and some analytical form near the origin 
must be assumed. In this work we provide then a method 
to extract information from the lattice QCD data, with- 
out assuming any special analytical form, as a dipole, 
tripole or exponential function. We start by requiring an 
overall consistency between the prediction of our model 
for the four form factors and the lattice QCD data for the 
physical n~ mass. Then, the model extends naturally 
down to the = limit, inferring the behavior of the 
il^ system in the region not covered by the lattice data. 
As a bonus relatively to an ad-hoc parameterization, we 
obtain even information about the microscopic structure 
such as the admixture percentage of each D- state and 
the momentum distribution of the wavefunction. 

The procedure assumes that meson cloud dressing is 
not significant in the D,^ system. Although the present 
lattice QCD results cannot rule out that sea quark dress- 
ing (meson cloud) may be important, they suggest that 
dependence in the light quark (pion mass) is small at 
least for the charge and magnetic dipole form factors. 
This topic will be investigated in the future. 

Our final results imply an unexpected large Dl state 
mixture (6.6%). The confirmation or disproval of this 
result will be possible once lattice data for Gm3 became 
available. More precise data for Gmi {Geo is already 
very precise), particularly at high Q^, where the weight 



of the D states is larger, will also be useful for an even 
better estimate of the both D-state admixtures. 

Our final result for the electric quadrupole moment of 
the Q- is Qa- = (0.96 ± 0.02) x IQ-^efm^. In the liter- 
ature the existing results for Qo- correspond to the in- 
terval (0.4-4.0) X 10-2efm2 [fl, [iMl III [S [II , and 
our result has a magnitude consistent with this range. 
The value extracted directly from lattice QCD data, as- 
suming an exponential dependence [s^l, corresponds to 
(1.18 ± 0.12) X lO-^efm^. Note that our result satisfies 
simultaneously the constraints of the three form factors 
{Geo, Gmi and Ge2), and therefore can be given with a 
smaller band of uncertainty. 

The numerical value of Qn- /sq- {gq- = —1), because 
it is a positive number, can be interpreted in a nonrel- 
ativistic formalism as the charge distribution being ex- 
tended and flattened along the equatorial region, as it 
was also predicted for the A+ baryon jlj. Note however 
that this interpretation has been questioned, and some 
authors suggest a different concept of deformation, based 
on the transverse electric quadrupole moment in the in- 
finite momentum frame [45144 7| . 

Comparing the fl~ with the A"*" baryon, in the covari- 
ant spectator quark model it is interesting to notice that 
the difference in the D3 admixture, 0.72% for the A and 
0.11% for the fi" does not correspond to a reduction of 
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G_B2(0) in the same proportion. This shows that the mo- 
mentum distribution in the overlap integral between S 
and D3 states is very different in both systems and it has 
to be taken into consideration. 

An intrinsic limitation of our calculation is the deter- 
mination of the form factors only in first order of the 
admixture coefhcients. The exact calculation of the form 
factors including the contributions from transitions be- 
tween D states is in progress Nevertheless, the ef- 
fect of these transitions can be estimated approximately 
by the correction implied by the deviation of the nor- 
malization constant N from 1. In the present case this 
correction is about 3.4%. 

Finally, in the future, a precise calculation of the 
magnetic octupole form factor of the n~ using lattice 
QCD will be very useful for a better understanding of 
the fl~ electromagnetic structure. The magnetic oc- 
tupole moment Oq- = GM3(0)2]gT, was estimated al- 
ready in Refs. [ill, [l3j E^- Our result corresponds to 
(1.27 ± 0.04) X 10~^efm3. Although the experimental 
determination of O^- may not be possible in practice, 
an evaluation based on QCD in a regime where meson 



excitations are expected to be small, may help to se- 
lect models. In particular for our model, it can allow a 
more accurate determination of the wavefunction struc- 
ture, which carry information on the shape of the electro- 
magnetic distributions inside that baryon, and also the 
the relative contribution of the D-wave components in 
the wavefunction of the . 
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